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ABSTRACT 

Unbounded outflows in the form of highly collimated jets and broad winds appear to 
be a ubiquitous feature of accreting black hole systems. The most powerful jets are 
thought to derive a significant fraction, if not the majority, of their power from the 
rotational energy of the black hole. Whatever the precise mechanism that causes them, 
these jets must therefore exert a braking torque on the black hole. Consequently, we 
expect jet production to play a significant role in limiting the maximum spin attainable 
by accreting black holes. We calculate the spin- up function — the rate of change of black 
hole spin normalized to the black hole mass and accretion rate — for an accreting black 
hole, accounting for this braking torque. We assume that the accretion flow onto a Kerr 
black hole is advection-dominated (ADAF) and construct easy-to-use analytic fits to 
describe the global structure of such flows based on the numerical solutions of Popham 
& Gammie (1998). We find that the predicted black hole spin-up function depends 
only on the black hole spin and dimensionless parameters describing the accretion 
flow. Using recent relativistic magnetohydrodynamical (MHD) numerical simulation 
results to calibrate the efficiency of angular momentum transfer in the flow, we find 
that an ADAF flow will spin a black hole up (or down) to an equilibrium value of about 
96% of the maximal spin value in the absence of jets. Combining our ADAF system 
with a simple model for jet power, we demonstrate that an equilibrium is reached at 
approximately 93% of the maximal spin value, as found in the numerical simulation 
studies of the spin-up of accreting black holes, at which point the spin-up of the hole 
by accreted material is balanced by the braking torque arising from jet production. 
The existence of equilibrium spin means that optically-dim AGNs that have grown 
via accretion from an advection-dominated flow will not be maximally rotating. It 
also offers a possible explanation for the tight correlation observed by Allen et al. 
(2006) between the Bondi accretion rate and jet power in nine, nearby, X-ray luminous 
giant elliptical galaxies. We suggest that the black holes in these galaxies must all be 
rotating close to their equilibrium value. Our model also yields a relationship between 
jet efficiency and black hole spin that is in surprisingly good agreement with that seen 
in the simulation studies, indicating that our simple model is a useful and convenient 
description of ADAF inflow - jet outflow about a spinning black hole for incorporation 
in semi-analytic modeling as well as cosmological numerical simulation studies focusing 
on the formation and evolution of galaxies, groups and clusters of galaxies. 

Key words: galaxies: jets; galaxies: nuclei; galaxies: active; accretion, accretion discs; 
black hole physics; magnetic fields 



1 INTRODUCTION 

Astrophysical black holes are characterized by just two prop- 
erties, their mass and angular momentum, since they will 
typically have zero charge (Blandford & Znajek 1977) . The 
angular momentum, J, of a Kerr black hole of mass M., can 
be defined in terms of a dimensionless spin j = Jc/GMl 
which must lie in the range —1 < j < 1. There is mounting 
evidence that astrophysical black holes are rotating and that 



the corresponding spin is a critical variable in establishing 
the observed properties of the accreting black hole systems 
(see discussion in McKinney & Gammie 2004; Gammie, 
Shapiro & McKinney 2004) as well as the magnitude and 
the impact that the outflows from these systems will have 
on their surroundings — see Nemmen et al. (2007) and ref- 
erences therein. The energy associated with this spin can 
be up to 30% of the gravitational mass in an extreme Kerr 
black hole (McKinney & Gammie 2004) and is, in princi- 
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pie, available for extraction (Penrose 1969) and so is a likely 
power source for phenomena such as jets from active galactic 
nuclei, microquasars and gamma ray bursts. Understanding 
how black holes acquire their spins and, in the cosmologi- 
cal context, what the distribution of spins is, is therefore of 
fundamental importance. 

Black holes can gain angular momentum either via 
mergers with other black holes or from the material that 
they accrete. In this paper, we restrict our discussion to ac- 
cretion. Bardeen (1970) was the first to consider the issue. 
Assuming that the accretion proceeded via a cold disk and 
that the angular momentum of the accreting material was 
always aligned with the spin of the black hole, he noted that 
an accreting black hole would spin up to its maximal rate 
in a finite time (i.e. after accreting around 1.5 times the ini- 
tial black hole mass). Thorne (1974) recognized that thin 
disks radiate and as a result, the black hole spin would be 
limited to j w 0.998 because the black hole would prefer- 
entially swallow negative angular momentum photons emit- 
ted by the accretion flow. A number of studies since then 
have suggested that this limit may not hold (or that the 
limiting spin may be higher) if accretion takes place via a 
geometrically thick flow. Popham & Gammie (1998), who 
numerically solved the equations describing the steady state 
radial structure of hot, geometrically thick, relativistic, op- 
tically thin, advection-dominated accretion flow, found that 
the equilibrium spin can be as high as 0.99997 or as low as 
0.7, depending on the efficiency of angular momentum trans- 
port in the flow; they did not, however, take into account ra- 
diation swallowing. In the case of pure advection flows where 
none of the energy dissipated in the flow is radiated away, 
this is not an issue. Finally, as noted by J. M. Bardeen (re- 
lated by Thorne 1974), the effects of magnetic fields could 
potentially limit the spin to a significantly lower value by 
exerting torques on the accreting material. 

The potential role of magnetic fields has received in- 
creasing attention in the past decade, beginning with the 
works of Krolik (1999) and Gammie (1999), who demon- 
strated that the magnetic fields expected on the basis of 
flux freezing should have order unity effects on the accre- 
tion flow and therefore, on the amount of energy released 
in that flow. The effects of magnetic fields have been stud- 
ied further with the general consensus that magnetic stresses 
can significantly alter the spin-up of the black hole (Li 2000; 
Li 2002; Agol & Krolik 2000). In particular, Agol & Krolik 
(2000) used the second law of black hole thermodynamics to 
place constraints on the maximum accretion efficiency for a 
given black hole spin. They found a maximum efficiency of 
0.36, which would be attained if the equilibrium spin was 
0.94. 

In recent years, the use of sophisticated general rela- 
tivistic, magnetohydrodynamic (MHD) numerical simula- 
tions (e.g., Koide et al. 2000; Koide 2003; de Villiers 
& Hawley 2003; McKinney & Gammie 2004; Gammie, 
Shapiro & McKinney 2004; Komissarov 2005; de Villiers et 
al. 2005; Hawley & Krolik 2006; Komissarov et al. 2007; 
Punsly 2007; Beckwith et al. 2008; McKinney & Bland- 
ford 2009) have yielded considerable insights into the role of 
magnetic fields in accretion flows, and especially the efficacy 
of magnetic stresses at extracting rotational energy from a 
black hole-accretion flow system. In the specific simulations 
presented by McKinney & Gammie (2004) and Gammie, 



Shapiro & McKinney (2004), stresses associated with mag- 
netic fields in the accretion flow resulted in the black hole 
evolving towards an equilibrium spin (i.e. dj/dt — 0) of 
j ~ 0.9 (specifically equilibrium occurs somewhere between 
their j = 0.90 and j = 0.94 simulations). Gammie, Shapiro 
& McKinney (2004) note that this equilibrium value is rel- 
evant only to geometrically thick accretion flows. Krolik, 
Hawley & Hiroise (2005) too reach a similar conclusion (i.e. 
that there is an equilibrium spin at approximately j = 0.9) 
by analyzing the MHD simulations of de Villiers & Hawley 
(2003). A more in-depth review of these studies can be found 
in Krolik (2005). To summarize, current MHD simulations 
of spinning black holes accreting from a geometrically thick 
accretion flow predict an equilibrium spin somewhere in the 
range j — 0.9-0.94, the exact value being uncertain due to 
numerical limitations and the small number of simulations 
available. At higher values of j, the black hole is spun-down. 

A key feature of all of these MHD simulations is the 
emergence of unbounded outflows. These outflows can occur 
in the form of a highly collimated, Poynting-flux dominated 
component generated both by the magnetized, rapidly rotat- 
ing accretion flow as well as the rotation of the black hole 
(de Villiers et al. 2005; Punsly 2007; McKinney & Bland- 
ford 2009), and as a broad, mildly relativistic, typically 
matter-dominated, component that originates in the accre- 
tion flow 1 . In the present paper, we will use the generic label 
"jet" to refer to the outflows, regardless of whether they be 
broad or narrow, matter or Poynting-flux dominated, and 
instead we will distinguish between outflow engendered by 
the magnetic fields anchored in the rotating flow (hereafter 
referred to as the "disk" jet) and that due to field lines an- 
chored on the black hole event horizon (the "black hole" 
jet). The black hole jet is primarily electromagnetic and the 
disk jet is hydromagnetic, where the energy (and the angu- 
lar momentum) is carried by an electromagnetic flux as well 
as a kinetic flux of matter. 

While a detailed, universally accepted, description of 
the mechanisms underlying the origin and power of these 
jets remains elusive, the combination of the results from sim- 
ulations with different parameters and analysis of simplified 
analytic models, suggests the following tentative picture that 
we use as a basis for our model for the jets (see Appendix 
B): 

As a starting point, we are interested in black holes 
situated at centres of galaxies and clusters of galaxies. In 
such systems, the gas accreting onto the black hole will 
have originally flowed in from kiloparsecs-scale via a cool- 
ing flow and a quasi-spherical accretion flow. During the 
latter stage, the magnetothermal instability is expected to 
organize the existing weak magnetic fields into primarily ra- 
dial fields (e.g. Sharma, Quataert & Stone 2008). We will 
therefore assume that there is a weak, large-scale poloidal 
component present in the accretion flow onto the black hole. 
Although the detailed geometric structure of the flow as it 
converges onto the black hole will depend on whether the 
gas is able to radiate away its thermal energy, generically 
the flow will be rotating and axisymmetric in character. In 



1 McKinney and collaborators refer to this component as the 
"disk wind" (c.f. McKinney 2006; McKinney & Blandford 2009) 
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the present paper, we will focus primarily on geometrically 
thick, advection-dominated accretion flows (ADAFs). 

Several MHD processes are expected to arise in the flow: 
the shearing of the poloidal field within the rotating fluid 
will give rise to a toroidal field and the flow will be subject 
to MHD turbulence. These features have been confirmed by 
the various simulations. The magnetic fields anchored in the 
rotating flow will give rise to a combination of outflow of 
material from the disk and its surroundings, and an elec- 
tromagnetic flux, with the former dominating if the black 
hole is spinning relatively slowly. As recognized by (Meier 
1999; see also Meier 2001), if the black hole is spinning, the 
frame-dragging of the inflowing gas within the ergosphere 
will enhance the magnitude of the outflow. In fact, accord- 
ing to Punsly (2007), not only do the numerical simulations 
by Hawley & Krolik (2006) confirm that the power of the 
hydromagnetic disk jet indeed grows with black hole spin, 
but that the rise in power is very steep at high black hole 
spin rates and that this steep increase is primarily due to 
growth in the power of the electromagnetic component of 
the disk jet. 

Finally, when the accretion flow reaches the black hole 
event horizon, the gas is expected to drain off the field lines. 
In the case of a non-rotating black hole, the magnetic pres- 
sure will cause the field lines, which for all intents and pur- 
poses can be thought of as being anchored on the black 
hole's event horizon, to establish a nearly radial configu- 
ration in the polar regions similar to the split-monopole 
structure first described by Blandford & Znajek (1977). 
In the event that the black hole is spinning, the winding 
of the magnetic field lines in the ergosphere will drive heli- 
cal twists along the magnetic tower that manifests as highly 
collimated, Poynting-flux dominated jets. 

In the picture outlined above, the combined disk+black 
hole jet therefore draws its energy from the gravitational 
energy released by the accretion flow as well as the rotational 
energy of the black hole itself. Consequently, the jet power 
does not vanish in the case of a Schwarzchild black hole. This 
is consistent with simulation results in that studies exploring 
the relationship between the outflow structure and power, 
and the spin of the black hole (e.g. McKinney & Gammie 
2004; de Villiers et al. 2005; Hawley & Krolik 2006) all 
find that disk winds persist even when the black hole is not 
rotating. However, in cases where the magnetic field lines are 
appropriately orientated 2 , the presence of a rapidly spinning 
black hole can greatly enhance the outflow power. This link 
between black hole spin and outflow power has long been 
indicated from theoretical considerations (c.f., Blandford & 
Znajek 1977; Punsly & Coroniti 1990; Meier 1999; Meier 
2001). 

Given that the jets gain a fraction of their power by 
tapping the rotational energy of the black hole, we can es- 
timate the braking torque exerted on the black hole and 
consequently, the maximum spin attainable by an accret- 
ing black hole driving powerful jets. Utilizing the under- 
standing of jet production developed by these simulations, 
we construct a simple, analytical model for the maximum 



2 Recent simulation studies suggest that an ordered polodial com- 
ponent is optimal (c.f., Beckwith et al. 2008; McKinney & Bland- 
ford 2009). 



spin attainable through accretion. To do this, we assume an 
advection-dominated accretion flow (ADAF) onto a spinning 
black hole whose global structure resembles the numerical 
solutions of Popham & Gammie (1998), and we examine 
how the accretion-driven spin-up of a black hole is modified 
when jets are produced as a consequence of that accretion. 
Using a simple model for the jet power, where a fraction 
of the jet power is extracted from the rotational energy of 
the black hole, we derive a modified spin-up function and 
thereby predict the equilibrium value of the black hole spin. 
In the process, we also determine the dependence of jet ef- 
ficiency as a function of black hole spin. 



2 SPIN-UP OF A BLACK HOLE 

We wish to compute the rate of spin-up for a black hole 
accreting from a geometrically thick, advection dominated 
accretion flow (ADAF) that also engenders powerful jets. We 
will begin by examining the classical calculation of the spin 
up of a black hole accreting from a thin disk and with no jet 
production. We will then proceed to modify this calculation 
to describe accretion from an ADAF and finally, quantify 
how jet production alters this spin-up. 



2.1 Spin-up by a thin accretion disk 

As stated above, we quantify the angular momentum of 
a Kerr black hole by the dimensionless parameter j (= 
Jc/GMl). We follow the notation of Shapiro (2005) and 
define a dimensionless spin-up function s(j) by 



s(j) = 



dj M. 

dt M. ' 



(1) 



where M._o is the rate of rest mass accretion. For a stan- 
dard, relativistic, Keplerian thin-disk accretion flow with no 
magnetic fields, we will denote the expected spin-up function 
due to accretion by so{j). Shapiro (2005) gives 



so(j) = Cisco — 2ji?isco, 



where 

-Cisco (j) = 
-Bisco(j) = 



Vnsco(r 2 sco - 2jyrisco + f) 
nsco \J r\ sco - 3nsco + 2j A /ri S co 



2nsco + jy^isco 



nsco \Afsco ~~ 3nsco + 2jVnsco ' 



(2) 

(3) 
(4) 



are the (dimensionless) specific angular momentum and spe- 
cific energy of the innermost stable circular orbit (ISCO) of 
the black hole respectively. This relationship implicitly as- 
sumes that (a) the inner edge of the disk coincides approx- 
imately with the radius of the innermost stable circular or- 
bit, inside of which the centrifugal force is unable to balance 
gravity and the gas begins to free-fall inward, and (b) that 
the torque at this inner boundary is negligible since the gas 
in the plunging region will quickly accelerate to supersonic 
speeds and lose causal contact with the material upstream. 
This "no-torque" boundary condition has been a subject of 
some debate but both Afshordi & Paczynski (2003) and 
Li (2003) have independently shown that in the case of 
thin disks, the assumption is reasonable. Moreover, recent 
simulations (Shafee et al. 2008) confirm that the torque is 
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negligible (at least for ] « 0). The radius of the ISCO orbit 
(in units of the gravitational radius GM./c 2 ) is 

nsco = 3 + A 2 (j) -y/[3- Ai{j)][3 + Mj) + 2A 2 (j)] (5) 
where 

Ai(j) = l + [(l-j a ) 1/3 ][(l+j) 1/3 + (l-j) 1/a ], 
Mi) = Vzf+Mj) 2 - (6) 

The net spin-up rate in the case of thin disk accretion 
is given by so(j) + s r ad, where s ra( j is the spin-down rate 
due to radiation swallowing as given by Thorne (1974) 3 . 
As illustrated in Figure 1, the net spin-up function for the 
standard thin-disk (green curve) is positive for all j < 0.998, 
and therefore lets the black hole spin up to j « 0.998 in finite 
time as noted by Shapiro (2005) 4 . 

2.2 Spin-up by an ADAF 

Thin disks are expected to be features of radiatively efficient 
accretion flows and therefore, typically associated with opti- 
cally luminous active galactic nuclei. The black hole systems 
that are of particular interest to us, i.e., supermassive black 
holes that reside at the centres of massive elliptical galax- 
ies, are typically optically-dim. In the inner regions, these 
systems are thought to be accreting in an optically thin, 
advection-dominated mode, in which the density of the ac- 
creting gas is sufficiently low that the gas cannot radiate effi- 
ciently; the energy dissipated during the accretion remains in 
the flow in the form of thermal energy and is eventually car- 
ried across the horizon (Esin, McClintock & Narayan 1997; 
Narayan, Mahadevan & Quataert 1998; Narayan, McClin- 
tock 2008). 

Whether the accretion flow is geometrically thin and 
radiatively efficient, or geometrically thick, optically thin, 
low luminosity and advection-dominated, or for that mat- 
ter, geometrically thick, optically thick, luminous but yet 
advection-dominated, depends largely on the mass accretion 
rate. Estimates of mass accretion rates onto the supermas- 
sive black holes in massive elliptical galaxies (c.f., Allen et al. 
2006) yield very low values and at low mass accretion rates, 
the flow is expected to be an optically thin, low-luminosity 
ADAF. 5 

The physical structure of an ADAF is very different 
from that of a thin disk. Consequently, the thin-disk expres- 
sion for so is unlikely to be correct for an ADAF. A priori, 
one would expect both the specific angular momentum and 
the specific energy of the accreted material to differ from 

3 Note that the quantity da„/dlnM plotted in Figure 6 of 
Thorne (1974) is defined using the actual change in mass of 
the black hole rather than mass flow rate in the accretion flow 
and therefore differs from our quantity s by a factor of -Eisco- 

4 Shapiro did not include the radiation breaking term and so 
finds that the black hole spins up to j = 1. We note that the 
radiation braking term is small for all j and discounting it would 
make no noticeable difference to any of the curves in this Figure. 

5 Since we are specifically interested in the properties of optically 
dim AGNs in massive elliptical galaxies, we will not consider the 
luminous, advection-dominated flows any further. Henceforth, we 
will use the label "ADAF" as a short-form for "optically-thin, 
low-luminosity, advection-dominated flows" . Such flows are also 
sometimes referred to as radiatively inefficient accretion flows. 



the thin-disk expectations. We therefore write the spin-up 
rate for ADAFs as: 

s' = £adaf0', a) — 2j_EADAF, (7) 

where £adaf and Sadaf are normally evaluated at the hori- 
zon. 

In the case of an ideal, unmagnetized ADAF, in which 
the gas starts out cold at large radii and all dissipated energy 
is advected along with the flow, the (dimensionless) specific 
energy of the accreted material is -Eadaf ~ 1 (Popham & 
Gammie 1998), while the specific angular momentum of the 
accreted gas is typically sub-Keplerian (Narayan, Kato & 
Honma 1997; Popham & Gammie 1998; Manmoto 2000). 
Recent simulations, which include the effects of magnetic 
fields and allow for the exploration of less constrained flow 
structures, confirm the expectation for the angular momen- 
tum, that it is typically less than the thin-disk result, but 
suggest that the energy of the accreted material is in fact 
close to the thin-disk value (see Figure 5 of McKinney & 
Gammie 2004 and also de Villiers et al. 2005). In this 
paper, we will show results for both values for the specific 
energy of the accreting material, the thin-disk value as in- 
dicated by the simulations as well as the ideal ADAF value. 

To quantify the accreted angular momentum (and, as 
described in the next section, to estimate jet power), we 
need a model for the ADAF accretion flow. To this end, we 
make use of the results of Popham & Gammie (1998), who 
numerically solved the equations describing the radial struc- 
ture of steady state, optically thin, advection-dominated ac- 
cretion flow in a Kerr metric (see also, Gammie & Popham 
1998). Though more cumbersome to use than the simple, 
easy-to-manipulate, self-similar solutions of Narayan & Yi 
(1995); Narayan, Mahadevan & Quataert (1998); Bu et. al 
(2008), they offer one key advantage over the latter: an im- 
proved treatment of the flow in the innermost regions of the 
ADAF where the self-similar results break down. The result- 
ing ADAF model includes the effects of magnetic fields to 
the extent that the flow is treated as a mixture of a frozen-in, 
isotropically tangled, turbulent magnetic fields and ionized 
plasma. It does not, however, account for large-scale mag- 
netic fields or MHD effects, like the amplification of mag- 
netic fields by dynamo-like processes and MHD winds. On 
the other hand, Popham & Gammie (1998) explicitly illus- 
trate how the flow structure depends on the black hole spin 
j, the adiabatic index 7 quantifying the strength of the tur- 
bulent magnetic fields in the flow, the viscosity parameter a 
governing the transport of angular momentum in the flow, 
and the advected fraction of the dissipated energy /. 

For the purposes at hand, we adopt their / = 1 (pure 
advection flow) results and to facilitate the use of these solu- 
tions in this (and future) work, we have derived analytic fits 
to key ADAF structural quantities, from which additional 
properties of the flow can be derived. These simple fitting 
functions are listed in Appendix A. For simplicity's sake, 
we shall assume that the rate of angular momentum accre- 
tion by the black hole per unit rest mass accreted, £adaf, 
is equal to the value of the specific angular momentum of 
the accreting fluid near the horizon. This is tantamount to 
assuming a "no-torque" boundary condition at the horizon. 
Strictly speaking, this condition does not apply to geomet- 
rically thick flows but in the model under consideration, the 
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Figure 1. Spin-up parameter of a black hole as a function of spin, j, for the case / = 1 and 7 = 1.44. Lines show results from the 
present work. The green line corresponds to spin-up by a thin disk (i.e. so) with the radiation swallowing term, s ra( j, included (but no 
braking term due to jet power). The blue line corresponds to spin-up by an ADAF (i.e. s' ) without the radiation swallowing term or 
jet braking term included, while the red line includes the jet braking term. We do not factor in the s ra( j term for the ADAF cases since 
these flows are radiatively inefficient and s ra[ j is small for all j (i.e. removing it would make no noticeable difference to the lines in this 
figure). The left plot uses a(j; E = 1) and the right plot is based on a(j; E = -Eisco)- The spin-up function for an ideal ADAF (E = 1) 
is shown as a solid curve while the results for E = -Eisco appear as dashed curves. Also shown are the simulation results from Krolik 
(2005) and Gammie, Shapiro & McKinney (2004) (blue and magenta points respectively). The horizontal line indicates the equilibrium 
state where s = 0. 



rate of angular momentum transport is very small near the 
horizon and the two quantities agree to within ~20%. 

In Popham and Gammie's treatment, the viscous stress 
is assumed to result from the combined effect of correlated, 
large-scale, time-averaged Reynolds and Maxwell stresses 
that they associate with MHD turbulence. They model this 
stress using a modified form of the a- viscosity prescription 
designed to ensure causal behaviour and their a parame- 
ter is defined in terms of the total pressure, which includes 
the contribution of the turbulent magnetic field pressure. 
Popham & Gammie (1998) assume that a, as well as all 
other parameters mentioned above, are constant over the 
radial extent of the flow. MHD simulation studies of magne- 
tized accretion flows suggest that this approach does not do 
justice to the rich dynamics involved. For one, MHD simu- 
lations of accretion flows around black holes have long indi- 
cated that the effective value of a is a strong function of ra- 
dius in the inner regions of the flow (Hawley & Krolik 2001; 
Hawley & Krolik 2002), rising from ~ 0.1 over the bulk of 
the flow to ~ 1 in the plunging region (see, for example, 
Fig. 4 in Hawley & Krolik 2002 and Fig. 12 in McKinney & 
Narayan 2007), with the large value of a corresponding to a 
non-negligible flux of angular momentum being carried out 
from inside the ISCO as a result of the presence of magnetic 
fields. In the present work, we recognize these limitations of 
the a-viscosity formalism and use a merely as a parameter 
characterizing the structure of the accretion flow. We choose 
its value so that the spin-up function that includes the ef- 
fects of outflows (see §2.3) is comparable to the results seen 
in MHD simulations. 



Allowing a to vary as a function of black hole spin j, 
we find (see Fig. 1): 

a(j) = 0.025 + 0.055j 2 for Badaf = 1, 

a(J) = 0.025 + 0.4/ for Badaf = -Eisco- (8) 

This dependence of a on j makes qualitative sense if one 
expects that the angular momentum transport is facilitated 
by magnetic fields and that the amplitude of these fields 
ought to grow as the black hole spin is increased. Such a 
correlation between magnetic field strength and black hole 
spin has been noted in numerical simulations (Hirose et al. 
2004). Interestingly, as we shall show below, the resulting 
flow structure, when combined with our jet power model, 
leads to jet efficiencies comparable to those seen in complex 
MHD simulations, suggesting that our model may be useful 
for calculating global properties of AGN systems that are 
required as input in theories of galaxy, group and cluster 
formation. We note that given our usage, our Q-parameter 
values cannot be directly compared to physical quantities 
such as the effective magnetic a of Hawley & Krolik (2002) 
and McKinney & Narayan (2007). In fact, the very concept 
of a- viscosity is not central to our argument. Any accretion 
flow with similar structural properties in the inner regions 
should yield similar results. 

The net spin-up function, s' , for our ADAF model is 
shown in Figure 1 as cyan curves. The solid curves show 
the spin-up function for an ideal ADAF (Eadaf = 1) 
while the dashed curves show the spin-up function for the 
case Eadaf = Eisco, as suggested by numerical simu- 
lations. The curves in the left plot are computed using 
a(j; Eadaf = 1) while those in the right plot make use 
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of a(j; -Badap = Sisco)- We note that in calculating the 
spin-up rate for an ADAF, we do not include a spin-down 
term due to radiation swallowing since radiation is not a 
factor in the case of pure advection flows. However, we com- 
puted the s ra d term out of curiosity and found it to be small 
for all j; including it makes no noticeable difference to our 
results. For the case -Eadaf = Ejsco (-Eadaf = 1), the 
spin-up curve is positive for j < 0.96 (0.97) and slowly spin- 
ning black holes will be spun up as accretion proceeds. For 
higher values of j, the spin- up curve is negative and the ac- 
cretion from the ADAF will spin down rapidly whirling black 
holes. The cross-over point between these two regimes de- 
marcates the equilibrium black hole spin value for accretion 
via an ADAF. Specifically, the equilibrium spin is j — 0.96 
for .Eadaf = -Eisco and j = 0.97 for Sadaf = 1. The lower 
equilibrium spin value relative to the thin disk case is the 
consequence of the accreting fluid's specific angular momen- 
tum being sub-Keplerian. 

2.3 Effects of jet production on spin- up 

We now consider what happens when a jet is launched in 
conjunction with accretion onto a black hole. There is a 
considerable body of work indicating that the launching of 
jets is most efficient when the accretion flow is advection- 
dominated (e.g. Meier 2001; Churazov et al. 2005) and 
that jet production is suppressed in thin disks (Livio et al. 
1999; Meier 2001; Maccarone et al. 2003). We will there- 
fore use the ADAF model described previously as our work- 
ing platform. Jets are, however, inherently magnetohydro- 
dynamic phenomena (Blandford & Znajek 1977; Blandford 
& Payne 1982; Punsly & Coroniti 1990; Meier 1999; Meier 
2001) and our ADAF model does not treat such effects. As 
a workaround, we follow the approach outlined by Nemmen 
et al. (2007) in developing a model for jet power and graft 
it onto our ADAF model. This approach implicitly assumes 
the basic structure of the accretion flow will not be signif- 
icantly altered by the implied presence of strongly ordered 
magnetic fields in the inner regions. Recent numerical sim- 
ulation results of Beckwith et al. (2008) seem to support 
this. 

The jet model that we use is described in Appendix B. 
The basic physical idea underlying our model is that the 
electromagnetic flux and plasma outflow that comprise the 
jet are due to a rotating helical tower of magnetic field lines 
engendered by the differential frame dragging of a preferen- 
tially poloidal field anchored on the event horizon, as well as 
the combined effect of differential rotation of the plasma in 
the body of the disk and differential frame-dragging of the 
plasma inside the black hole's ergosphere on the poloidal 
field lines anchored in the accretion flow (see Hirose et al. 
2004; Beckwith et al. 2008). In this model, the jet therefore 
draws its energy from the energy released by accretion as 
well as the rotational energy of the black hole itself, with 
the contribution from the latter dominating when the black 
hole is spinning close to its maximal rate. 

The fact that jets launched from the black hole- 
circumnuclear ADAF disk system can derive some or all 
of their power from the rotational energy of the black hole 
suggests that they will exert a braking torque on the black 
hole. Using the irreducible mass of the black hole (Misner, 
Thorne & Wheeler 1973; eqn. 33.58): 



~ 0.4 




Figure 2. The equilibrium spin, j cq (a) = j(s nc t = 0;a), as a 
function of a for the Eadaf = 1 case is shown by the thick solid 
curve while that for -Eadaf = ^ISCO case is shown as thin solid 
curve. The dots show the location of the equilibrium spins. The 
corresponding maximum radiative efficiency of the black hole, 
e rad[jcq(a)], is shown by the dashed lines while the jet efficiency, 
r\ = Pj c t/ M,fiC 2 is shown by the dotted lines. 



M > r 

M., irr = — 



(1+ 



1/2 



(9) 



we can derive how the spin j changes due to jet braking. 
From eqn. (9) we find 



4 / M., irT y 

. M.\ M. ) 



dM. Mm \ Mm J J 

The spin-down rate is then 

d? _ dj [/BH,«Pdisk,jct + -PBH.jct] 

dt ~ dM. c 2 

and the corresponding jet spin-down parameter is 
Mm dj 



(10) 



(11) 



Sjet 



M, o dt 



(i + v/i 



-ff+j 2 



/BH,.fdisk,jot + fBH.jct 



(12) 



In the above equations, Pdiskjet and PBHjet correspond to 
the jet power associated with the unbounded outflows en- 
gendered by magnetic fields anchored in the rotating flow 
and the black hole event horizon, respectively, and /bh,» is 
the fraction of the disk jet power that is extracted from the 
black hole's rotational energy. 

The disk and BH jet power are given by (see Ap- 
pendix B for details) 



Pdisk 



jet — 



v 2 r 



V 
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x(^ + ^) 2 > d3) 

and 

P BH , Jot = -(l-VgrM.v — ^^^T^—) 

x(» + ^) 2 , (14) 

where all quantities are evaluated at the ISCO for the disk 
jet and at the static limit for the BH jet, the relativistic 
boost factors j r and 7^ are defined by eqns. (A16) and (A17) 
respectively, (3^ = ^Jl — I/7J, the metric factors A and T> 
are defined by eqns. (A15) and (A14) respectively, V is the 
radial velocity of the flow given by eqn. (A2), £adaf is given 
by eqn. (A26), T is a dimensionless temperature given by 
eqn.(A19), and g is a magnetic field enhancement factor due 
to shearing (see Appendix B) and is given by g = exp(ii)r), 
where to is the angular velocity of spacetime rotation given 
by eqn. (A18) and r is a characteristic timescale available 
for field enhancement given by eqn. (B7) . 

The BH jet power vanishes as the black hole spin 
approaches zero. In other words, PBHjct = for the 
Schwarzchild black hole. However, the unbounded flows in 
the numerical simulations (de Villiers & Hawley 2003; de 
Villiers et al. 2005; Hawley & Krolik 2006) do not vanish for 
this case, indicating that the accretion flow itself also makes 
a contribution. The disk jet power does not go to zero as 
the black hole spin goes to zero, but it is not independent of 
the black hole spin either. It grows with increasing spin. The 
effects of frame-dragging of the accretion flow within the er- 
gosphere of a spinning black hole enhances the power of the 
disk jet by a factor of g 2 (see Appendix B). We use this to 
estimate the fraction of the power that must come from the 
accretion power. Since g = 1 in absence of frame-dragging, 
we adopt /bh,« = 1 — g 2 as a reasonable approximation 
for the fraction of the disk jet power that is extracted from 
the black hole's rotational energy. This approximation states 
that regardless of the value of the black hole spin, accretion 
will always provide some amount of power to the outflows. 
For a Schwarzchild black hole, the disk (and total) jet power 
is then driven entirely by accretion power (Mc 2 ) while for 
a maximally spinning Kerr black hole only a fraction of a 
percent of the disk and total jet power is extracted from ac- 
cretion power, the rest being extracted from the black hole 
spin. Specifically, 72.8% of the disk jet power (about 93% 
of the total jet power) is drawn from the black hole's rota- 
tional energy for j = 0.8, rising to 92.9% (96.6%) at j = 0.9. 
We note that the results that we present in this paper are 
not sensitive to the details of this approximation providing 
that, for j close to unity, the majority of the jet power is 
drawn from the black hole's rotational energy. As we will 
show below, this is indeed the case. 

The net spin-up function for Eadaf = 1 and Eadaf = 
Eisco ADAFs with jet outflows are plotted in Fig. 1 as red 
curves. Our results display the same qualitative behaviour 
as the Krolik (2005) and Gammie, Shapiro & McKinney 
(2004) simulation results, which we also show. This corre- 
spondence with the simulation results holds even if the a 
parameter is treated as a constant. Our results, therefore, 
make a strong prediction that there should be an equilib- 
rium value of black hole spin where we may expect to see 



the spin distribution be truncated (if black holes spin up 
from j = via an ADAF). 

For a(j) as described in eqn. (8), the spin- up func- 
tions for ADAFs with and without jets are nearly identi- 
cal at low spin values (j <C0.8). As noted previously, the jet 
power at low spin values is largely drawn from the accre- 
tion flow and is relatively low. At high spin values, however, 
the spin-up rate is noticeably lower once the jet power is 
taken into account, with the effect being more pronounced 
for the ideal ADAF case. Note that for j < 0.92, s nct is 
positive and low spinning black holes will be spun up as 
accretion proceeds, even if the accretion is accompanied by 
outflows. For j > 0.94, however, s nct is negative; the combi- 
nation of jets and angular momentum flux from the ADAF 
will spin down rapidly whirling black holes. The cross-over 
point between these two regimes demarcates the equilibrium 
black hole spin value for accretion via an ADAF. Specifi- 
cally, the equilibrium spin is j = 0.93 for Eadaf = 1 and 
j — 0.92 for Eadaf = Eisco The a value corresponding to 
the equilibrium spin value is 0.073 if -Eadaf = 1 and 0.317 
if -Eadaf = Eisco- 

The equilibrium spin depends on the dimensionless pa- 
rameters a, 7 and / that describe the structure of the 
ADAF. Setting aside the fact that we have chosen our a(j) 
so that the resulting spin-up function agrees with simula- 
tion results, one can ask how sensitive the equilibrium value 
is to the value of a at the point where s not = 0. In Fig. 2 
we plot the equilibrium spin j(s n et = 0) as a function of a 
(with fixed / = 1 and 7 = 1.444). As the plot illustrates, a 
decrease in the value of a results in a higher value of Jeqi 
and vice versa; the dependence, however, is relatively weak. 
Over the range, a — 0.05-0.5, the change in the equilibrium 
spin value is of order 10%. 

Fig. 2 also shows the run of jet efficiency, defined as 
T) — Pj e t / M.,oc 2 corresponding to the equilibrium spin value, 
where, as noted previously, M m ,o is the rate of rest mass 
accreting onto the black hole 6 . For equilibrium spin value 
j ~ 0.92, the jet efficiency factor is r\ = 0.06 and 0.16 for 
Eadaf = Eisco and Eadaf = 1, respectively. For values of 
a ;>0.3, the jet efficiency parameter is nearly constant. This 
is because the corresponding equilibrium spin values are not 
so high and the jet power comes largely from accretion. The 
equilibrium spin value simply represents the maximum spin 
state that accretion can drive the system into. The decline 
in equilibrium spin (and consequently jet efficiencies) with 
increasing a is due to the reduction in the angular momen- 
tum of the ADAF material. For a <J0.2, jet efficiency rises 
steeply with decreasing a. Even though the magnitude of 
the equilibrium spin value has not increased by much, the 
system is now moving to a state where the jet power is com- 
ing to be dominated by the rotational energy of the black 
hole. The equilibrium spin values seen in numerical simula- 
tions suggest that those systems are close to this latter state. 
In this state, the angular momentum being carried by the 
jet is not negligible and the equilibrium spin value is truly 
a result of a competition between spin-up due to accretion 
and spin-down due to the jet. 



6 The jet efficiency, r), defined in this way can exceed unity since 
the jet can draw power not only from the accretion flow, but also 
from the black hole spin. 



8 Andrew J. Benson & Arif Babul 



1 F 



0.1 - 



0.01 



0.001 - 



o.ooo: 



This work; E ADAF = 1 
This work; E ADAF = E ISC0 

^BH + ? 7disk 

^BH 

^disk 



/ - 

/ 

/ 

/ 



-7 



/ 



- 



Hawley & Krolik (2006) 
de Villiers et al. (2005; 







0.2 



0.4 



0.6 



0.1 



Figure 3. The jet efficiency, r) = Pj et /M, t oc , predicted by our ADAF model with 7 = 1.444 is shown as a function of j (lines). The 
solid line shows the results for our -Eadaf = 1 model while the dashed line shows the results for our Eadaf = ^ISCO- Magenta lines 
show the black hole (electromagnetic) jet efficiency, cyan lines show the disk (electromagnetic and matter) jet efficiency, and blue lines 
show the sum of the two. Estimates of ?7tot from the numerical simulations of Hawley &: Krolik (2006) and de Villiers et al. (2005) are 
shown by green and red circles respectively, with the filled points corresponding to the efficiency factor of the total jet power and the 
open circles corresponding to the efficiency factor of the Poynting-flux component. 



For completeness, we also show in Fig. 2, the maxi- 
mum radiative efficiency of the black hole, e ra d(snct = 0) = 
1 — -Eisco[j(snot = 0)] as a function of a. For both of our 
fiducial models, we find e ra d(snet = 0) « 0.16 and like the 
equilibrium spin, the magnitude of e ra d is only weakly de- 
pendent on a. We would be remiss if we did not mention 
that this radiative efficiency does not, of course, apply to 
an optically thin ADAF. However, we would expect this ef- 
ficiency to be relevant to an accretion system that has been 
spun up by an ADAF and then transitions to a radiatively 
efficient thin-disk mode but has not yet been in this new 
state long enough to alter the black hole spin. For complete- 
ness, we note that in general, the radiative efficiency of an 
AGN can be modeled approximately as 



, (m/m cr ) , if m < m c 

EADAP ~ CradU) X < . 

1, if m > m c 



(15) 



where rh is the Eddington-scaled mass accretion rate and 



m cr = 0.03 is the rate above which the accretion flow be- 
comes radiatively efficient (Merloni 2008). 

In Fig. 3, we show the jet efficiency as a function of j 
for our fiducial Sadaf = 1 and -Eadaf = Eisco (solid and 
dashed blue lines), and compare these to the available re- 
sults from numerical simulations. The magenta curves show 
the black hole jet efficiency, which vanishes at zero black 
hole spin; the cyan curves show the disk jet efficiency, which 
asymptotes to tyiisk ~ 0.002 as j vanishes; and the blue 
curves show the efficiency factor for the total jet power. The 
simulation results are shown as green and red circles. The 
filled circles show the efficiency of the combined power in 
electromagnetic and matter 7 outflows, and the open circles 
show the efficiency associated with the electromagnetic jet. 
While there remain significant differences between the var- 
ious numerical calculations of r\ (e.g. there is an order of 



7 Like both Hawley & Krolik (2006) and de Villiers et al. (2005), 
we remove the rest mass contribution to the matter outflow. 
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magnitude difference in the jet power at j — 0.5 as deter- 
mined by the two different calculations) , there is a clear and 
rapid increase in jet efficiency with rising j. This is repro- 
duced by our model, which is in reasonably good agreement 
with the simulation results. Specifically, the simulation re- 
sults indicate that even at j = 0, there are outflows from 
accretion flow-black hole system. Our model reproduces this 
non-zero jet efficiency at zero spin and indicates that it is 
powered entirely by the rotation of the accreting plasma. At 
intermediate spins (j « 0.5 — 0.8), the black hole jet domi- 
nates. At high spins (j > 0.9), the simulation results show 
a rapid rise in power. Our black hole jet model, which is 
based on the Blandford-Znajek model (Blandford & Zna- 
jek 1977), does not rise as rapidly and this has been one 
of the many criticisms of this particular model. Our total 
jet power, however, does show a steep rise and the results 
are consistent with the simulation results. The steep rise 
is driven by a sharp increase in the disk jet power. In our 
model, we cannot distinguish whether this steep rise in the 
disk jet power is due to increase in the kinetic flux of the 
matter outflow, the electromagnetic flux, or both. However, 
an analysis of the Hawley & Krolik (2006) numerical sim- 
ulations has led Punsly (2007) to conclude that for high 
spins, the effects of frame-dragging can not only result in 
disk jet that is Poynting-flux dominated but also enhance 
its power so that it dwarfs the black hole jet. 



3 DISCUSSION 

Jets appear to be a ubiquitous phenomenon in accreting 
black hole systems, with both the accretion flow as well as 
the black hole's rotational energy contributing to the jet 
power. The question of how much energy can be extracted 
from an accretion flow-rotating black hole system has been 
addressed previously by a number of groups, the most recent 
being Nemmen et al. (2007). The jets, however, also carry 
away angular momentum from the system and most of these 
studies have not considered the implications of this angular 
momentum extraction. Since jets emitted by systems with 
rotating black holes derive a fraction of their power from 
the spin of the accreting black hole (Meier 2001; Nemmen 
et al. 2007), they must exert a braking torque on the black 
hole. Jets, therefore, must play a central role in limiting the 
spin-up of black holes. 

In a nutshell, the spin equilibrium is the result of the 
coupling between the black hole, the accretion disk and the 
jet outflows mediated by magnetic fields. The detailed mech- 
anism through which the magnetic fields couple the rota- 
tional energy of the black hole to the outflows is a topic of 
a number of on-going analytic and numerical investigations, 
and we presented a summary of some of the recent insights 
and thinking about the origin of the jet phenomenon in §1. 
In this paper, we have described a simple calculation that 
attempts to capture some of the aspects of the complex pro- 
cess summarized above in the form of a highly simplified 
model. The resulting calculation, based on the ADAF model 
summarized in Appendix A and jet model summarized in 
Appendix B, indicates that any black hole which has under- 
gone significant accretion (i.e. more than doubling its mass) 
from an ADAF-type flow will have a spin of approximately 
j « 0.92 providing it remains undisturbed afterwards (i.e. no 



mergers or subsequent accretion by other means). Moreover, 
the predicted equilibrium spin should be independent of the 
initial spin state. For example, a black hole that begins ac- 
creting from an ADAF with an initial j > j(s nct = 0) (where 
the jet power is very large) will be rapidly spun-down to the 
equilibrium value. Essentially, any black hole driving a suf- 
ficiently powerful jet will reach an equilibrium spin j < 1. 
This fundamental result — that black holes driving powerful 
jets must experience a braking torque which plays a central 
role in limiting their spin — is independent of the details of 
the model that we employ. 

The calculation we present is in excellent agreement 
with results from MHD simulations of accreting black hole 
systems. Specifically, we show that by using a model that 
matches the spin-up function measured in the numerical 
MHD simulations, we are able to self-consistently explain 
the jet efficiencies measured in the same simulations. This 
dual agreement is a pleasant surprise, given the simplifica- 
tions and approximations inherent in our quasi-analytic ap- 
proach and suggests that our model represents a convenient 
useful description of ADAF inflow - jet outflow that can be 
easily incorporated in semi-analytic modeling as well as nu- 
merical simulations focusing on the formation and evolution 
of galaxies, groups and clusters of galaxies. Numerous stud- 
ies have shown that AGN feedback, both kinetic and radia- 
tive, is critical for preventing the formation of overluminous, 
much too blue, massive galaxies (Hopkins et al. 2005; Cre- 
tan et al. 2006; Bower et al. 2006; Somerville et al. 2008) 
and for tempering massive cooling flows that ought to have 
been present in at least 30% of galaxy clusters (McCarthy 
et al. 2008; Bildfell et al. 2008), but which the observations 
suggest are much weaker (Kaastra et al. 2001; Peterson et 
al. 2001). The nominal explanation is that radiative energy 
loss that underlies the cooling flow phenomena is being par- 
tially offset by heating by outflows from supermassive AGNs 
in cluster galaxies. 

Recently, Allen et al. (2006) analyzed X-ray and optical 
data from set of nine nearby, X-ray luminous giant elliptical 
galaxies that show evidence of jet outflow-related activity 
and found a remarkably tight correlation between the esti- 
mated Bondi accretion rate 8 and the estimated jet powers 
in these systems of the form: Pjot = r/BondiMBondiC 2 , where 
7?Bondi ~ 0.02. The results are consistent with those expected 
from the ADAF-jet model of the kind presented in the paper 
especially if one allows for the fact that not all of the mass 
flow at the Bondi radius will actually accrete onto the black 
hole. For example, if M #i o = /Meondi where / < 1, then 
the jet efficiency factor, as defined and used in this paper, 
is r\ — r/Bondi//- A careful analysis of the constraints on r\ 
imposed by the trend found by Allen et al. (2006) led Nem- 
men et al. (2007) to conclude that the resulting black holes 

Given a distribution of gas about a central black hole, the most 
simple configuration describing the accretion of the gas onto the 
the black hole is the Bondi flow model (Bondi 1952), which 
assumes a non-luminous central source and a spherically sym- 
metric flow with negligible angular momentum. The resulting 
Bondi accretion rate can be written as MB on di = Tt\c s pr A , where 
Tj\ = 2GM §i o/Cj is the accretion radius, G is the gravitational 
constant, M #i o is the black hole mass, c B is the sound speed of 
the gas at ta, p is the density of gas at Tj\ and A is a numerical 
coefficient that depends on the adiabatic index of the gas. 
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must be spinning rapidly (i.e. j >0.8). An examination of 
Fig. 3 shows that the use of our model leads to very similar 
conclusions. More importantly, though, Fig. 3 shows that jet 
efficiency is a strong function of the black hole spin. Over 
the range, 0.8 < j < f .0, n) varies by approximately two or- 
ders of magnitude. This suggests that if the spin of the black 
holes in massive elliptical galaxies is unconstrained and any 
value between 0.8 and 1 is equally likely, then one would 
expect to find Pj ot varying by as much as two orders of mag- 
nitude for a given value of Meondi and therefore, the Allen 
et al. (2006) plot should have resembled a scatterplot. The 
fact that it does not indicates that the black hole spins are 
not unconstrained but rather, span a tight distribution. We 
argue, based on the findings presented in this paper, that 
the most natural value for the black hole spins to cluster 
about is the equilibrium spin value of j £s 0.92. 

Constraining the black hole spins to j « 0.92, with the 
corresponding jet efficiency of r\ ~ 0.06, further implies that 
~ 30% of the material Bondi flow accretes on the black 
hole. This surprisingly large fraction indicates that the black 
hole is very efficient at capturing material within its Bondi 
radius, which for the systems in question is approximately 
10-20 parsecs. 

Finally, we emphasize that the calculation arguing for 
the existence of a spin equilibrium at j ~ 0.92 is predicated 
on the accretion disk-black hole system driving powerful out- 
flows from an ADAF. Black holes accreting via thin accre- 
tion disks will behave very differently. The jet power pro- 
duced by a thin disk is around three to four orders of mag- 
nitude smaller than that from an ADAF for j w 1 (Meier 
2001). The maximum |sj e t| = 4.1 x 10~ 5 for a standard thin 
disk occurs at j = 0.94. This is much less than |s ra d| which 
is approximately 0.016 at the same j. The jet will have little 
impact on the black hole spin and our model predicts that 
a black hole spun up by accretion from a thin disk should 
reach the limiting spin of j = 0.998 found by Thorne (1974). 
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APPENDIX A: ADAF FLOW ONTO A KERR BLACK HOLE 

In this Appendix, we describe our model for an advection-dominated accretion flow onto a Kerr black hole. Our model is 
based on the results of Popham & Gammie (1998), who numerically solved the equations describing the radial structure 
of steady state, optically thin, advection-dominated accretion flow in a Kerr metric (see also, Gammie & Popham 1998; 
Manmoto 2000). The accreting fluid is assumed to be a mixture of a frozen-in, isotropically tangled, turbulent magnetic 
field and ionized plasma, and the key parameters that determine the structure of the ADAF are (i) the dimensionless black 
hole spin parameter j, (ii) the adiabatic index 7 quantifying the strength of the magnetic fields in the flow, (iii) the viscosity 
parameter a governing the transport of angular momentum in the flow, and (iv) the advected fraction of the dissipated energy 

We make a number of simplifying approximations and assumptions. The first of these is that we will assume, like Gammie 

6 Popham (1998), that the angular momentum of the accreting fluid is aligned with the angular momentum of the black 
hole. Additionally, we will also focus our attention specifically on the structure of the flow close to the equatorial plane. 
Also, as noted in the text, we shall set / = 1; that is, we shall only consider pure advection flows. We have discussed our 
treatment of the viscosity parameter a in the text. The adiabatic index of the fluid will depend on the relative contributions 
of the thermal and magnetic energy densities to the total internal energy density, U to t = Knag + U gaB , of the fluid. If the 
magnetic field contributes energy density Knag = B 2 /8ir to the admixture, then the corresponding isotropic magnetic pressure 
is Pmag = Knag/3 (Narayan, Mahadevan & Quataert 1998). Assuming that the magnetic pressure contributes a constant 
fraction of the total pressure, P mag = (1 — /3)P to t, it is straightforward to show that (c.f., Esin 1997): 

8-3/3 / \ 1 \ 

7 =6^- (A1) 

Since < /3 < 1, 7 is constrained to range within [4/3, 5/3]. For a mixture in which the magnetic field is in equipartition 
(i.e. I7g as = Knag), P = 2/3 and 7 = 1.5. (As a point of clarification, we note that all thermodynamic quantities, including 
magnetic field pressure and energy density, are measured in the fluid's local rest frame.) 

In eqn. (Al), we have assumed that the ionized plasma can be described adequately as an ideal, classical gas. Specifically, 
we ignore any relativistic corrections that may arise at high temperatures. Finally, in the discussion below, we shall, for 
convenience, work in geometric units (G = c = M. — 1); the units of mass, length and time are M., GM./c 2 , and GM./c 3 , 
respectively. We also emphasize that Popham & Gammie (1998) (see Gammie & Popham 1998 for details) do not take into 
account the effects of any large-scale magnetic fields that may be present, nor do they treat MHD effects, like the amplification 
of magnetic fields by dynamo-like processes. Consequently, neither does the model presented below. 

In the reference frame of an observer who is corotating with the fluid about the black hole at the Boyer-Lindquist radial 
coordinate r, the radial velocity of the accreting material (c.f. Popham & Gammie 1998; Figs. 1-4) is well described by fitting 
function V(r) where 



V(r) = - 


y/l-(l- 2/reff + (j/reff) 2 ), 




where 






r c « = 


r n + *(r)(r - r h ), 




*(r) = 


V1V2V3V4V5, 




Vl = 


91og(9z), 




V2 = 


exp(-0.66[l - 2a e fi] log[a cff /0.1] log 




V3 = 


1 - exp[-z {0.160' - 1) + 0.76}], 




VA = 


1.4 + 0.29065(j - 0.5) 4 - 0.87560 - 


-0.5) 2 ) + ( 


V5 = 


2.3exp[40(j - l)]exp[-15nsco(z - 


Zh)] + 1. 


CKeff = 


a [l + 6.450(7 - 1.444) + 1.355(7 ~ 


1.444) 2 ] . 



(A2) 

(A3) 
(A4) 
(A5) 
(A6) 
(A7) 

!3j + 0.45035)[1 - exp(-(z - z h ))], (A8) 

(A9) 
(MO) 

In this series of equations, z = r/risco and z h — r h /risco, where r h (j) = 1 + \f \ — j 2 is the radius of the event horizon and 
HscoO) is the radius of the ISCO (see eqn. 5). 

Popham & Gammie's steady state ADAF is characterized by a constant rest-mass accretion rate, M. The combination 
of this accretion rate and the above velocity can be used, via the continuity equation, to compute the radial distribution of 
the rest-mass density, p(r): 



-M 1-V 2 

'< r > = i^wV— ' (A11) 

where H(r) is the characteristic angular scale of the flow about the equator, which according to Popham & Gammie (1998) 
is given by 

n 2 = (A12) 
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where T is the dimensionless temperature defined in eqn. (A19). In these equations, r\ = (P to t + f/tot + p)/p is the relativistic 
enthalpy of the accreting fluid, 

vl = ^ {f + [1 - (juj) 2 }£ 2 ADAF - [(j^f/AKjrVv) 2 - y r Vv[2£ A v AF uj^j 2 }/VA} , (A13) 
is the effective vertical frequency, 

V(r) = 1 - 2/r + (j/r) 2 (A14) 
and 

A = l+j/r 2 + 2j 2 /r 3 (A15) 
are relativistic (metric-related) factors, 



and 



*=V 1+ % (A17) 

are relativistic boost factors associated with the radial and tangential motions of the ADAF fluid with respect to the Boyer- 
Lindquist coordinates (Popham & Gammie 1998), 

^-i^i = ^L, (A18) 

is the angular velocity, in the same coordinate system, corresponding to the local spacetime rotation (frame-dragging) enforced 
by the spinning black hole, £adaf is the angular momentum of the accreting fluid, and Ptot is the fluid's total (thermal plus 
magnetic) pressure. 

Popham & Gammie (1998) do not give the total pressure explicitly. Rather, they plot the generalized dimensionless tem- 
perature of the accreting fluid, T(r, j, a, 7) = Ptot / p, which can then be used to compute the total pressure. This dimensionless 
temperature is well described by the fit: 

where 

ti(7) = -0.270278 7 + 1.36027, (A20) 

£2(7) = -0.94 + 4.4744(7 - 1.444) - 5.1402(7 - 1.444) 2 , (A21) 

t 3 (a) = 0.841og 10 a + 0.919 - 0.643 exp(-0.209/a), (A22) 

U(j) = [0.6365nsco(j)- 0.4828] [1.0 + 11.9 exp {-0.838n 4 soo (j)}]. ( A23 ) 

t 5 (j) = 1.444exp[-1.01rJ , S co(i)]+0.1. (A24) 
The relativistic enthalpy can also be expressed in terms of the dimensionless temperature: 

J ,wl + _X-T(r,j > a,7). (A25) 
7-1 

Here 7 is the adiabatic index of the accreting fluid. 

The one remaining quantity that appears in the above equations but has yet to be defined £adaf, the specific angular 
momentum of the accreting fluid. We determine this from the product j|£adap that Popham & Gammie (1998) plot (c.f. their 
Figs. 1-4). Our fit for this quantity is: 

77-Cadaf = (^adaf) 2 + [(^adaf)i + 10 ( " £adaf)6 ][1.15 - 0.03(log 10 a + 3) 2 37 ], (A26) 
where 

(^adaf)i = 0.0871nsco - 0.10282, (A27) 

(i)£adaf) 2 = 0.5 - 7.7983(7 - 1.333) 1 ' 26 , (A28) 

(j-;£adaf)3 = 0.153(r ISCO - 0.6) 30 + 0.105, (A29) 

(>)£adaf)4 = (j?£adaf)3(0.97 - 0.2996)(1.202 - 0.08[log 10 a + 2.5] 2 ' 6 ), (A30) 

(t?/:adaf)5 = -I.87 + 4.299 - 0.018 + 0.018(log 10 a + 2) 3 ' 571 , (A31) 

(r?/:ADAF)6 = (77rADAF)4{[(0.14(log 10 r) (??£ADAF)5 + 0.23)/(t7/:adaf ) 4 ] 10 + 1 } ' . (A32) 
These relations cumulatively completely specify the basic structure of a "simple" ADAF. 
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APPENDIX B: JET POWER FROM AN ADAF 

In this Appendix, we describe a model for computing the jet power arising from our ADAF, including the jet launched from 
the disk and the jet launched from the black hole. We will assume that these two jets arise from the same physical mechanism, 
differing primarily in the radii from which they originate. There is a considerable body of work indicating that the launching 
of jets is most efficient when the accretion flow is geometrically thick, advection-dominated flow (e.g., Meier 2001; Churazov 
et al. 2005) and least efficient when it proceeds via a geometrically thin accretion disk normally associated with radiatively 
efficient AGNs (Livio et al. 1999; Meier 2001; Maccarone et al. 2003). While a detailed understanding of the extragalactic 
AGN jet phenomena remains elusive, the combination of physically insightful analytic studies by a number of authors over 
the past three decades and recent sophisticated general relativistic, magnetohydrodynamic (MHD) numerical simulations are 
beginning to yield important insights. There is now a general consensus that jets are fundamentally MHD events. Since our 
ADAF model described in Appendix A does not treat MHD effects, we will develop a separate model for jet power that can 
then be coupled to our ADAF model. 

The jet model that we adopt follows the construction outlined by Nemmen et al. (2007), which itself is an improved 
version of the scheme first proposed by Meier (2001). We recognize that this disjointed approach is not fully self-consistent, 
but we accept this and other associated limitations in favour of a simple, easy-to-use, model. That our model jet results agree 
reasonably well with those found in MHD simulations gives us a measure of confidence. As in Appendix A, the description 
below will make use of geometric units where G = c = M. = 1, and we will, as before, seek to keep the model simple by 
adopting as input ADAF properties computed in the equatorial plane of the flow. 

As our starting point, we assume that both the disk and the black hole jet powers are given by (Meier 2001): 

Pjot = h P^.-M^M] 2 - ( B1 ) 

where all quantities are evaluated at a characteristic radius associated with the region within which the jet forms. For the jet 
launched from the black hole, we will assume a characteristic radius equal to the static limit (i.e. the edge of the ergosphere) 
in the disk plane, 6 — ir/2, i.e. 7? ata tic = 1 + y/l — p cos 2 6 = 2 in gravitational units. For the disk jet, we identify the 
characteristic radius with -Risco under the assumption that the jet originates from the region interior to the inner edge of 
the ADAF and that this inner edge coincides with f?isco- We acknowledge that there is an ambiguity associated with this 
identification. As discussed by Krolik & Hawley (2002) and Watari & Mineshige (2003) , there is no reason for the inner edge 
to occur precisely at the ISCO and, furthermore, the location of the "inner edge" will depend on the physical property that 
defines that edge. The various simulation studies suggest that the launch region could be a factor of 2-3 inside the ISCO. 
Under our scheme, a smaller characteristic radius would lead to slightly higher jet power. For example, at j — 0.9, the jet 
power is enhanced by a factor of two if the innermost edge occurs half-way between the ISCO and the event horizon. We 
consider such order unity uncertainties as an acceptable compromise. 

In eqn. (Bl), Bpoi.oo i s the large-scale poloidal magnetic field, as seen by an observer at infinity in the Boyer-Lindquist 
coordinate system and fioo is the angular velocity of the rotating magnetic fields (and the plasma), as seen by the same 
observer. For the disk jet, this angular velocity takes into account the fact that the fluid itself is rotating as well as the effects 
of frame-dragging if the black hole itself is spinning; the black hole jet depends only on the angular velocity associated with 
the frame-dragging enforced by the spinning black hole: 

o ( r \-[ + w ( r ) for disk J et ' 

Moo(r) - j w(r) for Wack ^ jet (Hi) 

where ui is the angular velocity associated with frame-dragging (see eqn. A18) and 



~ £ ADAF (r) 1 m± {m) 
r 2 lr(r)i4,(r) \ A A (r) 

is angular frequency of the rotating fluid with respect to the local inertial observer, also often referred to as the zero angular 
momentum observer (ZAMO). This is an observer who is dragged in azimuth by the spinning black hole and orbits with 
angular frequency w. Gammie & Popham (1998) refer to this observer's reference frame as the locally non-rotating frame. 

To complete our calculation of the jet power, we need to estimate the magnitude of the large-scale poloidal component of 
the magnetic field, as seen by the observer at infinity in the Boyer-Lindquist coordinate system. This field is related to that 
in the fluid rest frame by a boost factor 7^: 

Bpol.oo = 7oc-Bpol, fluid (B4) 

where (Takahashi, Rilett, Fukumura & Tsuruta 2002): 

7oo = 7>-70 

(2jf3 4> /r 2 + VV)/VA. (B5) 

In our ADAF model, we assumed that the magnetic field is isotropically tangled. The field in an accreting flow, however, 
is unlikely to be isotropic due to the presence of weak, ordered magnetic fields, and the results of MHD dynamo effects 
induced by turbulence as well as large-scale shearing motions. Livio et al. (1999) (see also references therein) suggest that 
under such circumstances, the relationship between the toroidal and the poloidal components of the magnetic field can be 
approximated as Bpoi,H u id ~ H-B^,, fluid- The resulting magnetic energy density is (7 m a g = (1 + H 2 )B^ fluid /87r. Identifying this 
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with the magnetic field energy density in our ADAF (i.e., !7 m ag = 3(1 — /3)Ptot) and noting that Ti ~ 1/2 at i?isco (Popham 
& Gammie 1998), we can estimate B^ auid . The approach described here is based on the implicit assumption that the basic 
structure of the accretion flow will not change significantly by the implied presence of strongly ordered magnetic fields in the 
inner regions. Recent numerical simulation results of Beckwith et al. (2008) suggest that this is a reasonable assumption. 

Additionally, Meier (1999) has argued that in the inner regions of the ADAF, the shear associated the differential 
dragging of reference frames will drive a dynamo that will draw on the poloidal component of the magnetic field to generate 
and amplify the toroidal component at the expense of the black hole's rotational energy. This effect too has been tentatively 
observed in recent MHD numerical simulations (Hawley & Krolik 2006) . Denoting the corresponding field enhancement factor 
as g — exp (lot), 

24tt , 



V-^po^oo / 



7o 



-(l-/3)PtotS 



(B6) 



The amplitude of this enhancement factor depends on the black hole spin through the angular velocity of spacetime 
rotation, uj. In the case of a non-rotating black hole, to = and g — 1 (i.e., no field enhancement). The enhancement factor also 
depends on the characteristic timescale, t, over which the field amplification would occur. The timescale typically associated 
with MHD processes is the orbital period of the rotating fluid as determined by the local inertial observer, r = ~ Q~ , but 
doing so involves an implicit assumption that that radial motion of the fluid is negligible. Such an approximation is reasonable 
if the characteristic radius of the jet is sufficiently far from the horizon; however, in the case of the black hole jet, the static 
limit is always close to the horizon. For the disk jet, as the black hole spin parameter tends towards unity, J?isco approaches 
the horizon, the radial velocity grows rapidly and the flow becomes increasingly radial in character. The rapid inflow will 
decrease the amount of time available for field enhancement. We take into account the resulting correction by modifying the 
characteristic timescale for field enhancement as follows: 



1,^ 



1 

n' 



s/vv 



(B7) 



In the above equation, Tf ~ (r/VV)/v r is the characteristic inflow timescale in the local inertial frame. 



